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I. INTRODUCTION 

The computation of expectation values in statistical 
mechanics and in quantum field theory in its functional 
integral formulation requires taking averages of functions 
with a large number of variables. In the continuous case 
this means integration over manifolds with large dimen- 
sions. In such cases standard numerical integration tech- 
niques are no longer efficient and one has to resort to 
Monte Carlo methods. Unfortunately, in some applica- 
tions of great practical interest such as lattice quantum 
chromodynamics in the presence of a baryonic chemical 
potential [1, 2], the Boltzmann weight to be used in the 
averages is not positive or even real. This fact prevents a 
straightforward application of the Monte Carlo method 
in these cases. 

For a (real and positive) probability density function 
(PDF) P{x), expectation values can be estimated by a 
pure Monte Carlo method in which the points are in- 
dependently sampled from P{x). The dispersion in the 
estimate of (/) p in the pure Monte Carlo method is 
'^p{f)/ \^ where crp(/) is the variance of / and N is 
the number of independent sampling points [3, 4]. If 
sampling P{x) is very costly it may be more convenient 
to use the so called reweighting method, in which an aux- 
ihary PDF Po{x) is sampled instead, making use of the 
identity 



(/>i 



(PPn')^ 



(1.1) 



For a generic observable f{x) this method is less efficient 
as it suffers from the importance sampling problem: less 
points fall in the relevant region (i.e., where P{x), rather 
than Po[x), is large), and the dispersion ap{f)/ ^ Ncs in- 
creases. A typical signal of this problem is the presence of 
large fluctuations in both the numerator and the denom- 
inator in (1.1) when Pi^{x) and P{x) are too different. 

Nevertheless, the reweighting method can be trans- 
lated immediately to the cases in which the weight func- 
tion P{x) is complex. (With some abuse of language 
such complex density functions are still referred to as 



PDFs.) A positive Pq{x) is chosen (a standard choice be- 
ing^ Pq{x) — A|P(a;)|) and used to generate the samples. 
Although the importance sampling problem is not well 
defined here in general (being P{x) complex, sampling 
it is meaningless), the same large fluctuation problem 
which appeared in the case of positive P{x) is present 
here as a sign (or rather phase) problem, since often the 
average of the phase, {P/\P\)\p\ , is very small but not its 
variance.^ 

There are cases, however, where a sampling of a com- 
plex P{x) can be given a meaning. A clear instance of 
this is a PDF defined on R such that P{x) = Pq{x — ia), 
where Pq (x) is positive (for real x) and analytic on a re- 
gion containing M and R — ia. Then for /(z) analytic on a 
region containing R and R -I- ia, {f{x)) p = {f{x + ia)) p„ . 
In this case the Monte Carlo method can be applied us- 
ing as sampling points Zk = Xk+ia, k — 1, . . . , iV, where 
the Xk are generated from Pa{x), and taking as observ- 
able f{z) (the analytic extension of f{x)). This example 
suggests a technique consisting of sampling the complex 
plane, or more generally a suitable complexified version 
of the original manifold where the complex P{x) is de- 
fined, using an appropriate real and positive p{z), and 
trade the computation of {f{x))p by that of {f{z))p. In 
practice this program has been applied by means of the 
so called complex Langevin method [12, 13]. For posi- 
tive P{x) this approach produces a random walk which 
asymptotically samples the probability density function. 
This property is easily shown from the associated Fokker- 
Planck equation of which P{x) is the stable stationary 
solution. The stochastic differential equation can be ap- 
plied to the complex case using P{z), the analytical ex- 
tension of P{x). For some P{x) it can be shown that 
the random walk will reach a stable equilibrium which 
samples a certain p{z) with the correct expectation val- 



^ Taking Poix) to be proportional to |P(a;)| needs not be the op- 
timal choice in practice, see e.g. [5]. 

^ See [1, 2] for alternative techniques, such as Taylor expansion, 
in the context of lattice QCD and [6—11] for other approaches to 
the complex action problem. 
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ues [14]. In those cases the method is very useful since 
then cssentiahy the same local algorithms used in the 
real case apply. Unfortunately, for many relevant com- 
plex P{x) there is no steady solution or, if there is, it 
does not display the correct expectation values [15-31]. 
As shown in [27] the latter problem arises from the fact 
that the Fokker-Planck equation admits more solutions 
in the space of distributions than in the space of ordi- 
nary functions. However, the complex Langevin method 
guarantees at most that the PDF obtained from projec- 
tion of p{z) on the real axis is a distributional solution of 
the Fokker-Planck equation, not that it should coincide 
with the original P{x). In essence the problem is that 
in that approach everything depends on P{z) and this 
function does not favor integration along the real axis 
from integration along many other curves on the com- 
plex plane (which can be viewed as a real axis after an 
analytic change of variables). Regrettably, unlike the real 
Langevin case, no practical criterion is known to decide 
a priori whether a given P{x) will produce a steady ran- 
dom walk (rather than a state that looks stationary on a 
finite time computer simulation) and if so, whether such 
steady solution is really sampling the input PDF P{x). 
(See however [32, 33].) 

Given the limited success of the complex Langevin 
method, a natural question is to what extent a given 
complex weight can be represented at all by means of 
any suitable ordinary (i.e., positive) PDF on the com- 
plexified manifold. A partial answer was given in [34] 
where it was shown that one dimensional complex PDFs 
admit such representation. In [35] it was shown that all 
complex weights on R'' of the type Gaussian times poly- 
nomial of any degree and in any number of dimensions 
are also representable by positive PDFs on C''. Since this 
set is dense in L'^{R.'^) this suggests that representability 
is a quite general property. In the present work we ex- 
tend those results by proving representability for smooth 
complex PDFs on W\ which are either rapidly decreasing 
at infinity or periodic. The latter case is then general- 
ized to smooth PDFs defined on any compact matrix Lie 
group manifold. 



II. REPRESENTATION FOR PDFS OF THE 
TYPE CHARGE-DIPOLE PAIR 

Let P{x) be the complex "probability" density func- 
tion for which we want to compute expectation values of 
"observables" f{x), 



positive measure dii{z), in such a way that 



(/)P := 



jdtx{x)P{x)f{x) 
J dn{x)P{x) 



(2.1) 



P and / are complex in general, d^j,{x) is a positive mea- 
sure and the normalization J d^{x)P{x) is different from 
zero. A positive representation of P{x) is a real nonnega- 
tive probability density function p{z) defined on the com- 
plex extension of the original real manifold, with some 



(/>P = := 



Jdp{z)p{z)f{z) 
Jdfi{z)p{z) 



(2.2) 



where /(z) stands for the analytical extension of f{x). 
Of course a precise definition requires to specify the space 
of allowed test functions /. For a real manifold such as 
M'^ a usual choice is that of the set of polynomials. For 
complex probabilities which are periodic (real manifold 
equivalent to a torus) the natural test functions are of the 
type e*" '". In general, the larger the set of test functions 
the smaller the number of representable complex proba- 
bilities. Hereafter we will assume that P{x) and p{z) are 
normalized, i.e., the denominators in (2.1) and (2.2) are 
unity. 

The key idea of the construction is to decompose P{x) 
as a sum of simpler PDFs 



P{x)=J2Pn{x) 



(2.3) 



all of them with positive normalization, or, equivalently, 

P{x)^J2'^r.Pn{x) (2.4) 
n 

with normalized P„ and Wn > 0, and then find positive 
and normalized representations Pn{z) for the Pnix) so 
that 



(2.5) 



provides the desired positive representation of P{x). 

A suitable choice is to take the P„ of the form Dirac 
delta distribution plus a derivative of a Dirac delta. If P 
is viewed as a hypothetical (complex) charge distribution 
with positive total charge, it can be decomposed as a 
sum of simpler distributions of the form positive charge 
(delta distribution) plus a dipole (derivative of delta) at 
the same point, with complex dipolar moment. In 
this is 



Q{x;h) = d{x) + h-V5{x) , xG 



he&. (2.6) 



As shown constructively in [35] such distribution admits a 
positive representation of the form Gaussian times poly- 
nomial. A similar but simpler solution is considered here. 
Let us start by considering the one- dimensional PDF 



Qi{x) = 5{x) + S'{x), xe 



(2.7) 



It can be represented on the complex plane by the posi- 
tive and normalized PDF 



qi{z) 



1 

Stt 



1 



zeC, (2.8) 



that is 



Ic 



d^zqi{z)z'' 



dxQi{x) a;" 



0,1, 2,... (2.9) 
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where (P z = dRc z dim z is the Lebcsguc measure on . 
This can be checked by direct computation using polar 
coordinates, or alternatively, writing q\ as 

g{z) := ^e-l^/2|% (2.10) 

integrating by parts and noting that {z"')g = Sn,o- 

Now a positive representation of Q{x; h) is easily ob- 
tained as follows: 

q{z;h)= I d^ziqi{zi)S{z-zih). (2.11) 
Jc 

{S{z — zq) being the Dirac delta distribution at zq of the 
measure cf'^'z.) Indeed, 

(/W>«. = / d'''zf{z)q{z;h) 

Jc 

= [ (fzfiz) ( d^ziqi{zi)6{z-zih) 

Jc Jc 

= / d'^zif{zih)qi{zi) 
Jc 

= [ dxif{xih){5{xi) + S'{xi)) 
Jr 

= fio)-h-vm 

= [ d^xf{x)Q{x;h) 
Jr-^ 

= {fix))Q.- (2.12) 

In the fourth equality we have used that f{zih) depends 
analytically on zi and qi{z) is a representation of Qi{x). 

The support of the distribution q{z: h) is at most a 
plane, and Gaussian-like on that plane. From the Monte 
Carlo point of view this is considerably more efficient 
than the solutions presented in [35] which had a larger 
size (Gaussian-like with non vanishing width in all 2d 
directions). 



III. PERIODIC PDFS 

Next let us apply this technique to find positive rep- 
resentations for periodic PDFs functions on M.'^. The 
advantage of the periodic case is that the manifold is 
effectively compact, topologically a d-dimensional torus. 
Specifically we take P{x) with period 1 in each direction, 

P{x + n) = P{x), xeW^, neZ"^, (3.1) 

and normalized to unity with measure rfAt(x) = d^ X on 

the (i-dimcnsional torus T'^ = [0, 1]'' (suitably compacti- 
fied). In addition, we assume P{x) to be a smooth (i.e. 



C°°) function on the torus. A positive representation of 
P{x) is a smooth nonnegative function p{z) on which 
is periodic (in the real directions) 

p{z + n)=p{z), ZG&, neZ'^, (3.2) 

and such that (e^""'^) p = (e2""-^)p for ah n S Z'^, that 
is 

/ d-^a; P(a;)e2™-^ = / d^'^zpiz)e'^''''^-'' , n G 
Jj'd Jfd 

where the integration manifold of z is the complexified 
torus T"' = [0, 1]'' X R''. This manifold is non-compact 
in the imaginary direction. For subsequent use, we note 
that the periodic version of Q{x;h), 

QP{x; h)= Qi^ + ^) (3-4) 

admits the positive representation 

qP{z;h)= J2 q{z + n;h), (3.5) 

the sum being convergent (as a distribution). All the 
functions considered, being periodic, are well defined 
(one- valued) on T'^ or f^. 
We can decompose P{x) in the form 

P{x) = Po{x) + V • F{x) (3.6) 

where Po{x) is chosen to be everywhere strictly positive, 
smooth (on the torus) and normalized to unity. A sen- 
sible choice is Po(x) = 1, however we will keep the pos- 
sibility of a more general choice to maintain the analogy 
with the non compact case, below. By Hodge's decompo- 
sition theorem [36], there is a complex vector field F{x) 
which is also smooth on the torus (smooth and periodic 
on M'*). Such F{x) is non unique. A natural solution is 
just the electric field-like one: let 

pix) := Pix) - Po{x) , (3.7) 

then (3.6) becomes 

V • F{x) = p{x) . (3.8) 
The electric field-like solution is 

F{x) = V^x) , (3.9) 
where ^{x) is a smooth periodic solution of 

V2$(a;) = p(x) . (3.10) 



Smoothness is invoked to simplify the treatment. Likely the con- 
structions presented here can be extended to suitable spa<;es of 
distributions. 
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More explicitly, 



27rin-a: 



$(x) = J2 



1 Pn 2-Kin-x 



F{x) = y JL^e2-^"-==. 



(3.11) 



Note that F{x) is smooth because p{x) is smooth (i.e., 
\Pn\ decreases for large n faster that any inverse power). 
Po vanishes due to / d'^xp{x) = 0, since P{x) and Po{x) 
are both normalized to unity. 

Having a valid F{x), (3.6) can be rewritten as 



P{x) = I d^yPoiy) {d{x -y)+ H{y) ■ V5{x - y)) 

(3.12) 



with 



H{y) 



Fjy) 
Po{y) ■ 



(3.13) 



The vector field H{x) is also smooth on T'^. In turn this 
can be rewritten as 

P{x) = [ dSPo{y)Q{x-y;H{y)) 

d''yPo{y)QP{x-y;H{y)) (3.14) 



where QP{x:H{y)) is just Ql{x) of (3.4) with h = 
H{y). Therefore, P{x) has been written as a sum of 
charge-dipole distributions of the type (2.6). A positive 
representation of P{x) is thus readily found: 

p{z) = [ d%Po{y)q''{z-y;H{y)) 

dSPo{y)q{z-y;H{y)) (3.15) 



Jw 



where q^{z;H{y)) or q{z;H{y)) are just qP{z;h) or 

q{z; h) with h = H{y). The positive representation p(z) 
is also smooth and periodic on and so smooth on T"^. 



IV. POSITIVE REPRESENTATIONS FOR PDFS 
ON COMPACT LIE GROUPS 

Let us now consider the extension of the previous con- 
struction to Lie group manifolds. Such manifolds ap- 
pear often in applications such as lattice gauge theories. 
Specifically, as real manifold we take a connected and 
compact Lie group G of dimension d, which for conve- 
nience will be taken as a matrix group. Hence all ele- 
ments are of the form g = exp(a; • T) where x e K"^ and 
the d matrices T define a basis of the d-dimensional Lie 
algebra. The complexified manifold is the connected but 



non-compact Lie group G = {g = cxp(2; • T), z G C*}. A 
compact Lie group G admits a two-sided invariant metric 
g^j^ which endows G with a Ricmannian manifold struc- 
ture [37]. The metric can be normalized so that the cor- 
responding volume element, d'^x^/g, is the Haar measure 
of G normalized to unity [38].^ This measure is used 
in the evaluation of expectation values. For expectation 
values on G we also take its Haar measure. 

A natural set of test functions, is that of arbitrary 
polynomials of the matrix elements of ^ G G, /(<?) = 
/({.9y })- The analytical extension of these polynomials 
corresponds just to replace g with g, i.e., /({.9y })• 

The construction of a positive representation p{g) of a 
smooth normalized complex PDF P{g) on G is as follows. 
Let Pa{g) be strictly positive, smooth and normalized (for 
instance Pq = 1), then the diS'erence P — Pq integrates 
to zero on G and, by Hodge decomposition in a com- 
pact manifold, we can write (using differential geometry 
notation [36]) 



P{g) = Po{g) - d^F{g), 



(4.1) 



where F{g) = Fn{g)dx^ is a one-form on G. Equiva- 
lently. 



P{g) = Po{9) + ^^F^{g), 



(4.2) 



where V^u is the covariant derivative on G as a Rieman- 
nian manifold, with the usual Levi-Civita connection, 
andF^(5) = gt'" {g)F^{g). The smooth vector field F^(c/) 
is not unique and the electric field-like solution can be 
adopted, for instance, [36]: 



F^{g) = W{VY\P{9)-Po{9))- 



(4.3) 



Using (4.2) and upon integration by parts, 

{m)p = I dti{g)P{g)f{g) 

= I dt,i9)Poig){l~H^ig)d^)f{g),{4A) 



with 



H^{g) = Po{g)-'F^{g) . 



(4.5) 



Clearly, iJ^ (5)9^/(17) describes an infinitesimal point 
transformation, and this can be implemented as an in- 
finitesimal left translation by means of a smooth field 
H((?) taking values on the Lie algebra: 



(4.6) 



g = dct g^^, and arc just any local coordinates, not necessarily 
the normal coordinates appearing in e"''^. 
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With the help of the distribution Qi(xi) in (2.7), the 
expression in (4.4) can then be rewritten as 

{m)p= [ dn{g)Poig) f cia;iQi(a;i)/(e-^«(«)ff)(4.7) 

J Jr 

Since the dependence on xi in / is analytic, we can use 
the positive representation of Qi{xi) on C, qi{zi), to 
write 

{f{9))p = J dt,{g)Po{g)J^d'zm{zi)f{e^''^^'^g) .(4.8) 

This aUows to express the expectation value using a pos- 
itive representation on G, 

{f{g))p = {f{9))p ■■= J dn{g)p{g) f{g) , (4.9) 

with 

p{g) = J dfi{g)Po{g)q{g,g;n{g)) (4.10) 

and 

'l{~9,9;n{g))= [ cfzm{zi)S{e^''^^^^g,g). (4.11) 

Here S{go,g) denotes the Dirac delta distribution at 
for the Haar measure diJ.{g) on G. The extension to non 
connected G is obvious. 

V. POSITIVE REPRESENTATIONS FOR PDFS 

ON R'^ 

For PDFs on R'^ the previous construction based on a 
superposition of charge-dipole pairs can also be carried 
out, but it is technically more involved due to the lack 
of compactness of the real manifold. We will assume the 
normalized complex PDF P{x) to be in Schwartz space, 
i.e. smooth (infinitely differentiable) and rapidly decreas- 
ing at infinity {P{x) and all its derivatives go to zero 
faster than any inverse power of |a;|). The space of test 
functions can then be chosen as the set of polynomials of 

X. 

We choose a suitable positive and normalized Po{x) 
and write P{x) as 

P{x) = Po{x) + V • F{x) . (5.1) 

Proceeding formally we then obtain a positive represen- 
tation with 

p{z)= [ d'^xPoix) [ d^ziqi{zi)S{z-x-ZiH{x)), 

jR'' JC 

(5.2) 

with H{x) = F{x)/Po{x). 

There is a number of issues to be considered in this 

construction, such as the existence of a suitable Po(x) 
and of the vector field F{x), the convergence of p{z) as 



defined in (5.2), since the integral is on a non compact 
manifold, and finally, the convergence of the momenta 
|z|" of p{z) for all non negative n. 

Assuming by the moment that suitable Po{x) and 
F{x) exist, the two latter points will be fulfilled a fortiori 
provided the set of integrals 

/„= / d^''zp{z)\z^\ (5.3) 
Jc 

exist, where 2;" := n^=i -^m" '^ii^i nonnegative 
integers. Equivalently, the integrals 

/„= / d'^xPoix) [ d^Ziqi{zi)\{x + ZiH{x)r\ 

(5.4) 

should be convergent. We will take Pq{x) also in 
Schwartz space and momentarily assume that H{x) is 
smooth. In this case the only problem of convergence 
may can from the large \x\ or large \z\\ sectors of the 
integral. Since Pq{x) and gi(zi) are rapidly decreasing, 
convergence is ensured provided H{x) is bounded by a 
polynomial. This condition turns out to be rather re- 
strictive since it implies that F{x) should be rapidly de- 
creasing and moreover it should go to zero at a rate not 
much slower than Po(x) itself. 

Since both P{x) and Po{x) are smooth and rapidly 
decreasing, so is their difi^erence, p 

p{x) := P{x) - Po{x) , (5.5) 

and the equation on F{x) is 

V • F{x) = p{x) . (5.6) 

The electric field-like solution, (3.9), exists and is smooth 

but unfortunately, it will not be rapidly decreasing at 
infinity, in general: Because p{x) is in Schwartz space 
so is its Fourier transform, p{k). Then, {k/\k\^)p{k) is 
rapidly decreasing but not necessarily smooth at fc = 0, 
correspondingly, its Fourier transform is smooth but not 
rapidly decreasing, in general. Indeed, unless p{x) has 
radial symmetry, the electric field will have the (in three 
dimensions) well-known multipolar contributions which 
fall as an inverse power for large x [39]. To find a suit- 
able F{x) it is better to view this quantity as (minus) a 
dipolar density so that its divergence is the charge den- 
sity p{x). Because the total charge carried by p{x) is 
zero it should be possible to build it as a superposition 
of dipoles with complex dipolar moment. Moreover, if 
p{x) is smooth or rapidly decreasing or of compact sup- 
port, it should be possible to choose F{x) with the same 
properties (the support of F{x) being larger than that 
of p{x) in general). A suitable solution is as follows 

d\\'^~^p{\{x- a) + a) , x^a 

(5.7) 

where the basepoint a is an arbitrary point in K''. In- 
tuitively this solution corresponds to an arrangement of 
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linear chains of dipolcs, all chains starting at a and end- 
ing at each of the points in the support of p, in which 
each "positive" charge of one dipole is canceled by the 
"negative" one of the next dipole in the chain. The total 
charge at a equals that in p{x) which is zero. 

It is immediate to verify that V • Fa{x) = p{x) for all 
X ^ a. Also it is clear that Fa{x) is rapidly decreasing 
at infinity (or of compact support if p{x) is) and smooth 
at all points except the basepoint. At x = a the solu- 
tion Fa{x) will not be smooth in general, however this 
problem is easily fixed by taking a smooth average over 
a, 

F{x) = - [ dW^-^f cfiaC{a){x-a)p{X{x-a)+a), 

(5.8) 

where we choose the weighting function C{a) to be 
smooth, of compact support, positive and normalized to 
unity. Upon taking the average over a, the Fourier trans- 
form F{k) can be shown to be rapidly decreasing, in ad- 
dition to smooth, so F{x) belongs to Schwartz space. ^ 

The remaining issue is whether Po{x) can be chosen 
so that H{x) = F{x)/ Po{x) is polynomially bounded. 
Since the discussion is rather technical, this is shown in 
the Appendix. Here we only note that Po{x) should not 
go to zero too quickly for large \x\, as compared to P{x). 
Otherwise the falloff of p{x) and F{x) would be domi- 
nated by that of P{x) and the ratio F{x)/Po{x) could 
fail to be bounded by a polynomial. 

VI. ALTERNATIVE CONSTRUCTIONS 

The positive representation of a given complex PDF is 
by no means unique [35]. For periodic density functions 
we present here alternative positive representations to 
those studied in section III. They have the virtue of 
being quite localized on or T'^. This is convenient for 
their use in Monte Carlo integration since the localization 
decreases the fluctuations in the average over samples. 

The construction is based on decomposing the peri- 
odic and smooth P{x) as a weighted sum (with positive 
weight) of simpler PDFs of the form 

Pnix) = 1 + Ane^'''^-^ , neZ<^, nj^O. (6.1) 
® Of course, here we are addressing the generic case. In many 
practical cases this complicated construction is not needed. For 
instance, if p{x) is of the type Gaussian times polynomial, it is 

I 



where 2/1,2/2 € M'' are only constrained to satisfy 

2 _ g27r-n. j/i _ g27rn-3/2 



If the normalization of P{x) is positive such decompo- 
sition is clearly always possible after a discrete Fourier 
decomposition. Pn{x) has expectation values 

{1 for m = 
An iovm = n , meZ^. (6.2) 
otherwise 

Actually the parameter A^ is redundant. Indeed, the 
expectation values of two PDFs p{z) andp'(z) = p{z+a) 
are related by 

^g-27rim.2^^, = Q-2^im-a^^-2nim-z^^ ^ ae&. (6.3) 

Therefore, it is sufficient to find positive representations 
for An = 1, since other values are generated by a trans- 
lation (barring the trivial case An = 0): 

P„(a;) = l-Fe^"'"-^, n^O. (6.4) 

Let pn{z), with z = x + iy,hc a. positive representation 
of Pn{x). Because Pn{z) is also periodic with respect 
to a;, it can be decomposed in discrete Fourier modes, 
g27rjfc a; rpj^^ expectation values of e~27rsr7i z indicate that 

the Fourier modes fc = and k = n should be present. 

The mode k = n must also be present, for Pn{z) to be 
real. The minimum required is thus 

Pn{z) = /io(y)+/i(2/)e'"^"-'"+/i*(2/)e-2-"-»=, z = x+iy, 

(6.5) 

and /io(y) real. Furthermore, 

ho{y) > 2\h{y)\ (6.6) 

ensures the positivity of Pn{z)- 

A solution is easily found using the ansatz 

h{y) = ai5{y - yi) - a26{y - y2) , ai,a2>0. (6.7) 

At least two delta distributions with weights of opposite 
sign are required to satisfy {e^'^^'^'^)p„ = 0. Saturation 
of the bound provides ho{y) = 2\h{y)\. The remaining 

conditions (l)p„ = (e^^'^*"'^)p„ = 1 give the solution 

immediate, by going to Fourier space, to find a valid F{x) which 
is also of the type Gaussian times polynomial. 



(6.8) 

I 

The support of Pniz) is R"^ + R'^ (or rather T'^ + T'^). 
This is quite localized as compared to or T"^. The lo- 



Pn{z) 



2e27rn-i/i cos2(7rn • x)6{y - yi) + 2e'^'''^ '"^ sin^(7rn • x)diy - 2/2) 

g27rn-i/i _j_ g27rn-j/2 
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calization within this family of PDFs increases (meaning 
lower Shannon entropy relative to T'^ -\- T"^^) by taking 
n • 2/2 ^ — CO so that the branch y = yi dominates the 
density function. The limit does not exist, even in the 
weak sense, since the test functions e"*'" are of rapid 
growth on away from the real axis. 



VII. CONCLUDING REMARKS 

We have shown that very general complex PDFs can 
be represented by positive representations upon analyt- 
ical extension of the original manifold. This allows to 
sample them as required in the Monte Carlo method. 
The existence of these positive representations was not 
granted given the repeated failure of algorithms such as 
the complex Langevin approach when applied to general 
complex weight functions. A virtue of our construction 
is that it does not depend on the analytical extension of 
P{x) itself. In general one cannot expect P{z) to be well 
behaved on the complex manifold. This problem affects 
severely the complex Langevin approach. Nevertheless, 
it should be clear that our construction, as it stands, is 
more of formal interest than of practical use. First, the 
normalization of P{x) is generally not known; the nor- 
malization is needed to choose Po{x) and so to obtain 
the difference p{x). Second, given p{x), F{x) is not easy 
to construct. In fact, in practical cases, such as large lat- 
tices, only algorithms that are local (i.e., not much more 
than nearest neighbors) have a chance to be viable. This 
does not seem to be the case of, say, (5.8). 

Another limitation is related to importance sampling 
which was one of the main problems appearing in the 
reweighting method, (1.1). Strictly speaking this prob- 
lem does not exist if a pure Monte Carlo method is ap- 
plied to p{z). However, a lack of importance sampling 
manifests itself as an enhancement in the fluctuations of 
averages over samples. Because the positive representa- 
tion p[z) is not unique there are different representations 
of a given P{x), all of them with the same expectation 
values on analytic functions (but of course with differ- 
ent expectation values for arbitrary, non analytic, test 
functions). For instance, the convolution of a given rep- 
resentation p{z) with a positive function C(z) with ra- 
dial symmetry and rapidly decreasing at infinity yields a 
new representation p'{z) (this is because analytic func- 
tions are invariant under such convolution) [35]. The 
new representation will be wider, i.e., less localized, than 
the original one. As a consequence, although the ex- 
pectation values of (analytic) observables will be equal, 
{f)p = {f)p', the dispersion will be different, being larger 
for p'{z). Therefore, it becomes crucial in the represen- 



Thc Shannon entropy of a probability density function P relative 
to another PDF Po is defined as (— \og{P / Pq)) p . The larger the 
entropy the less localized is the probability. 



tation approach to find representations as localized, with 
entropy as small, as possible. 

This can be seen in another way. Eq. (5.2) indicates 
that the Monte Carlo method can be applied in the fol- 
lowing manner. First, a sample of x is generated from 
Po{x), and then, for each x the field H{x) is computed 
and a sample of zi is generated from qi{zi) to compute 
{f{x + zxH{x)))q^. This cc-dependent average is then 
itself averaged over the sample of Pq{x) to finally yield 
{f{z))p. However, in most, if not all, cases it will be more 
sensible to compute the average over qi exactly, 

{fix + z^Hix))),, = fix) - Hix) ■ Vfix) , (7.1) 

rather than using Monte Carlo. So a better approach is 
to directly compute 

{f)p = {f-H-Vf)p,, (7.2) 

(Interestingly, this formula does not require the analyti- 
cal extension of fix).) Such an approach is to be com- 
pared with the standard one, (1.1), 

{f)p = {f + Po-\P-Po)f)po (7.3) 

where, as in (7.2), we have assumed Pq to be normalized 
(not the case in practice). The two constructions are 
rather similar (formally H = P(^^V~^(P — Po)) so a 
priori there is no compelling reason to expect that the 
importance sampling problem gets better in (7.2) than 
in (7.3). [An exception would bo perhaps the cases in 
which / is particularly flat, since then this good property 
is enhanced by the derivative in (7.2)]. 

Despite these critical remarks, it remains the fact that, 
as shown here, the representation problem admits a solu- 
tion for very general complex PDFs P{x). Clearly, rep- 
resentability is a sine qua non condition for the success of 
any other approach based on analytical extension, what- 
ever the method used in the construction of the positive 
representation. 

From the mathematical point of view, several interest- 
ing problems pose themselves. One is finding more gen- 
eral construction methods, in addition to the one pre- 
sented here based on the chargc-dipole decomposition. 
Another is carrying out the construction for more gen- 
eral manifolds and more general complex PDFs. Also 
challenging is the problem of finding positive represen- 
tations of minimal entropy for a given P{x): as men- 
tioned before, by convolution it is always possible to in- 
crease the entropy, but there is no general mechanism 
to decrease it, so a minimum value is to be expected. 
[A similar entropy minimization problem has been found 
in a different context in [40].] If the positivity condi- 
tion on p{z) is relaxed, a quite localized representation 
is p{z) = Pix)5{y). Quite likely, this is the optimum so- 
lution if Pix) is positive. In the general case of complex 
P(a;), imposing piz) to be positive will probably imply a 
greater delocalization on the complex manifold. Finally, 
one can try to extend the representativity problem. In 
fact in our discussion we have made use of a map K from 
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the space of test functions f{x) to that of functions f{z) 
by means of analytical extension, as well as the adjoint 
map which is a projection from the set of PDFs p{z) 
on the complex manifold to that of complex density P{x): 
K\f) = I/), (iCtp| = (P| (note that the map p ^ P \s 
single- valued although P ^ pis not). That is 

(/)P = = {K^p\f) = {p\Kf) = {p\f) = if), 

(7.4) 

From this point of view, more general representations K 
could be sought which could be of practical interest. 
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APPENDIX A: POLYNOMIAL GROWTH OF 

H{x) 

Here we want to show that Pq (x) can be chosen in such 
a way that H{x) is bounded by a polynomial. We take 

it as evident that for any complex Schwartz function P' 
on K'^ there is a positive function Pq such that (i) Pq is in 
Schwartz space, (ii) Pq{x) > \P'{x)\ for all x e M'', and 
(iii) Pq is a decreasing function, i.e., Po{xi) < Po{x2) 
if > \x2\-'^ We apply this property to the function 
P'{x) := \x\'^'^P{x) and take Po{x) \x\-^'^P^{x) for 
\x\ > i?, for a sufficiently large R. The definition of Po{x) 
is completed in the region \x\ < R so that it is smooth 

To obtain Pq consider the auxiliary function 

P,','(cc)=sup{|P'(y)|,|y|>|£c|}. 

By construction Pq' satisfies all conditions (i-iii), except perhaps 
the property of being smooth. It seems obvious that Pq' can be 



and normalized. Therefore, for |a;| > R, \P{x)\ < Po{x) 
and 

|p(a:)| <2Po(a:;), 

Po{y)<^r^Po{x), \y\>\x\. (Ai) 



\y 



(As in Section V, p — P — Pq.) In order to use these 
inequalities we recall that in (5.8) the support of the pos- 
itive function C(a) is compact, so a lies inside a ball of 
some size Rc, and moreover A > 1. Thus the inequalities 
(Al) apply at the point X{x — a) + a provided R > Rc- 
This allows to write 



|F(a;)| < / dXX''-^ [ d'^aC{a 



Xttt^ . .^. Po{x), |a;|>P.(A2) 



\X{x - a) + a\^'i 
On the other hand |A(a; — a) + a\ > X{\x\ — \a\), hence 

/•OO P 

\F{x)\ < 2 dXX-^-^ / d'^aC{a) 



\x - a\\x\'^'^ 



< 



2 {\x\+Rc)\x 
d {\x\-Rcf'^ 

As a consequence, 



2d 



-Po{x), |a;|>P.(A3) 



\F{x)\ < k\x\Po{x), for \x\>R, 
2 1 + Rc/R 



k = 



d{l- Rc/Rf'^ 



(A4) 



and H{x) is bounded by k\x\. 

suitable smoothed out to obtain a valid Pq. 
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